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Abstract 

This conference proceeding presents the first results of the full CMS Tracker alignment based on several 
million reconstructed tracks from the cosmic data taken during the commissioning runs with the detector in 
its final position and magnetic field present. The all-silicon design of the CMS Tracker poses new challenges 
in aligning a complex system with 15 148 silicon strip and 1440 silicon pixel modules. For optimal track- 
parameter resolution, the position and orientation of its modules need to be determined with a precision 
of about one micrometer. The modules, well illuminated by cosmic ray particles, were aligned using two 
track-based alignment algorithms in sequence in combination with survey measurements. The resolution in 
all five track parameters is controlled with data-driven validation of the track parameter measurements near 
the interaction region, and tested against prediction with detailed detector simulation. An outlook for the 
expected tracking performance with the first proton collisions is given. 



1. Introduction 

Silicon tracking detectors in general purpose de- 
tectors like the Compact Muon Solenoid (CMS) at 
CERN [T] are built to reconstruct charged parti- 
cles trajectories (tracks). In a magnetic field, they 
are described by a helix. The track parameters 
are the curvature 1/pt (expressed as inverse trans- 
verse momentum), the impact parameters dxy and 
dz in the xy plane and along the principal axis of 
the experiment respectively and the polar angles 9 
and Their precise and accurate determination 
are paramount for the operation of tracking detec- 
tors with spatial resolution of the order of 10 /xm. 
Therefore the position of the modules needs to be 
known to better than this precision, which can be 
achieved by high mounting precision, stable frame, 
survey measurements and track based alignment. 
This article describes the track based alignment of 
the CMS inner tracker and the results obtained us- 
ing cosmic ray particles. Brief statements will be 
made on the use of survey information. 
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-"^The CMS coordinate system is defined as follows^: The 
origin is at the nominal collision point, the x-axis pointing 
to the center of the LHC, the y-axis pointing up and the 
2:-axis along the anticlockwise beam direction. 9 is measured 
from the positive z-axis and from the positive x axis. The 
radius r denotes the distance from the z-axis. 



1.1. The alignment problem 

Track based alignment can be described as a least 
squares minimization problem where the data from 
hits generated by tracks are used. A single resid- 
ual Tij for hit i along track j is the three dimen- 
sional distance between the predicted hit location 
from the track model and the physical hit informa- 
tion from the modules, calculated using the current 
knowledge of the geometry. Together with the co- 
variance matrix V the expression to be minimized 
is given in equation ([T]): 

tracks hits 

X^(P,q)= ^^'^%iP,(i])^r^Y,j{^,Cij) (1) 

3 i 

where p denotes the alignment parameters describ- 
ing the current geometry and denotes the track 
parameters of the j^^ track. In principle, this can be 
solved using standard techniques like solving nor- 
mal equations. 

The inner tracker at CMS consists of 1440 sili- 
con pixel modules and 15 148 silicon strip modules 
(figure [T]). Each module has six degrees of freedom, 
described in local coordinates u,v,w with respect 
to the geometric center of the module and rotations 
a, /3, 7 around these axes. In total we have to de- 
termine 99 528 parameters. For a typical alignment 
of the CMS inner tracker, around 10^ to 10^ tracks 
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are required, depending on which hierarchy levels 
(modules or larger units) are selected as objects to 
be aligned. Therefore the number of parameters to 
be determined in this procedure becomes at least 
of the order 0(10''). Solving it within hours, as re- 
quired for prompt alignment, is beyond the limit 
of the capabilities of the computers available to the 
experiment. 

2. Alignment algorithms used 

Two alignment algorithms were used to produce 
the results reported later in this article. Both aim 
to reduce the complexity of the problem so that 
it can be solved within hours on standard CPU'fEl 
They are distinguished by their scope: 

2.1. Global algorithm 

Solving the full alignment problem would pro- 
duce estimates for the alignment parameters and 
the track parameters. In our case, only the first are 
of interest. Restricting the solution to the align- 
ment parameters reduces the complexity to O(IO^) 
in our case. Using a clever scheme for setting up 
the matrix of the normal equations, this can be 
achieved using block matrix operations while us- 
ing the full information from the track parameters. 
This is implemented in Millepede-II[Z , an algo- 
rithm widely used for alignment purposes. Its ad- 
vantages are that it takes all correlations between 
modules and higher hierarchies into account. The 
algorithm works in a single step, however due to 
rejection of tracks with too high (outlier rejec- 
tion) , a small number of iterations are still required. 
The implementation in the CMS software frame- 
work uses a simplified helix model, in which ma- 
terial effects due to dE/dx are taken into account, 
but multiple scattering is currently ignored. This 
is a major disadvantage, as it limits the maximum 
resolution obtainable. At the time of this study, a 
memory limit allowed for an alignment of 46 340 pa- 
rameters at maximum in one step. An incremental 
procedure to align parts of the detector was used to 



^AU calculations were carried out on a batch farm 
at CERN consisting of nodes having 2 KSi2k on av- 
erage (KSi2k: Standard Performance Evaluation Cor- 
poration benchmark of Kilo Specmarks Integer year 
2000, http://wMW.spec.org/cpu2000/). Job parallelization 
among computer nodes is used whenever suitable and rea- 
sonable. 



overcome this limitation. To give an idea of the per- 
formance, typical time consumption (start-to-end) 
for a full alignment was about 4 hours. 

2.2. Local algorithm 

By assuming no track parameter dependence - 
dropping correlations between alignment parame- 
ters between modules - the problem can be reduced 
to solving the equation for single modules. Corre- 
lations between modules are recovered by iteration. 
The residuals r^- are calculated as the distance be- 
tween the physical hit data and the impact point 
from the track using the reconstruction procedure 
without the hit in consideration. This is imple- 
mented in the HlP-algorithm^. The major advan- 
tage of the implementation is the use of the same 
track model as in the track reconstruction (Kalman 
filter) and therefore all material effects are taken 
into account. On the other hand this algorithm 
convergences very slowly when the start geometry 
is not sufficiently close to reality. Typical time con- 
sumption for a full alignment was about 5 hours. 

2.3. Combined operation 

Both algorithms make use of job parallelization 
on the computer cluster for data collection steps 
(typically up to 100 computer nodes, no intercom- 
munication among concurrent jobs), while final cal- 
culations are carried out on a single machine. As 
the approaches are complementary, we used a com- 
bined method to benefit from the strength of both 
algorithms and to overcome their weaknesses. 

1. The global algorithm started from design ge- 
ometry. This resolved global movements and 
ended up in a geometry close enough to reality 
for efficient operation of the local algorithm. 
Despite the fact that the global algorithm is 
capable of aligning on several hierarchical lev- 
els simultaneously, the already mentioned pa- 
rameter limitation required the splitting into 
several steps. 

2. The local algorithm started from the outcome 
of the global one and resulted in a refined ge- 
ometry. 

Some of the plots in the result section will show the 
outcome of the individual algorithms together with 
the combined approach. A detailed description of 
all steps involved can be found in 2J. 
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Figure 1: Schematic view of the CMS inner tracker. The tracker consists of several subdetectors. The innermost part is 
the pixel detector (a barrel and two cndcaps at each side) surrounded by two barrel shaped strip detectors (TIB: tracker inner 
barrel, TOB: tracker outer barrel) and the endcap structures (TID: tracker inner disks, TEC: tracker endcap). [l] 



2.4- Survey information 

Survey data may come from optical surveys and 
coordinate measuring machines and are usually col- 
lected prior to or during installation. Alignment 
constants from such operations can be used as 

1. starting points for the alignment. This may 
enhance the convergence of an alignment al- 
gorithm, but the survey information looses its 
weight after the very first iteration. 

2. additional data for the alignment algorithm. In 
the local algorithm, this can easily be done by 
extending the sum of equation ([T]). The resid- 
uals for that are calculated as the difference 
between the position from the survey and the 
current reference geometry. 

Only the local algorithm used survey information 
in the results presented here. 



3. Results from commissioning with cosmic 
rays 

The results presented here are based on data col- 
lected in autumn 2008 during a period of cosmic 
ray data taking with a magnetic field of 3.8 T in 
the tracker volume. The total number of events 
detected by CMS during this campaign was about 
300 million, of which 3.2 million have hits in the 
tracker suitable for alignment. A cut on px > 
4GeV/c has been applied. The rate was about 



5 Hz. The fraction of tracks crossing the pixel detec- 
tor was 3% in the barrel and 1.5% in the endcaps. 
Data used for alignment and validation were not 
statistically independent due to the limited num- 
ber of events collected. Several low- and high-level 
approaches have been used to estimate and validate 
the alignment performance. 

CMS is designed primarily for tracks originat- 
ing from the nominal intersection point (including 
tracks from displaced vertices) and not for cosmic 
rays. For the tracker, this means that alignment 
is limited to parts with sufficient illumination from 
cosmic particles. We are also prone to deformation 
modes of the tracker which leaves the invariant. 
A known case is an elongation of the tracker along 
the z-axis, which is difficult to align using cosmic 
tracks only. 

3.1. Track distribution 

For each track of a data sample, the track 
is calculated. This is merely the second sum in 
equation ([ij , weighted by the number of degrees of 
freedom (ndof). 



■^track 

ndof 



1 



hits 

— ^rf(p,q)V-ir.(p,q) (2) 



A histogram of the distribution of these Xtrack 
lows for a low-level evaluation of the alignment. 
The results are shown in figure [2j where the im- 
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Figure 2: Track x^/ndof distribution. This plot shows 
the distributions for the unaUgned (dotted) tracker and after 
ahgning using the combined approach (soUd). For compari- 
son, the results after the alignment using the local (dashed- 
dotted) and global methods (dashed) are included. 



provement from the unaligned to the ahgned detec- 
tor is clearly visible. The combined approach shows 
the best alignment performance. 
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3.2. Distribution of the mean of the residuals 
(DMR) 

For each hit in a track of a data sample, the 
residual is calculated between the predicted posi- 
tion from the track and the actual hit where the 
hit has been removed from the track reconstruction 
in order to be unbiased by the hit under consid- 
eration. Such distributions were obtained for all 
modules individually. These are dominated by two 
effects: (1) track extrapolation uncertainties due to 
multiple scattering and (2) hit position uncertain- 
ties coming from the hit reconstruction algorithms. 
Both effects being random, they average out to a 
good approximation to zero if data from a suffi- 
cient number of hits is available. Misalignment is 
a systematic effect on these distributions. There- 
fore we determine the median (the 0.5 quantile of 
a distribution) for each module's distribution in or- 
der to measure such a systematic bias. These are 
then histogrammed for each subdetector, restricting 
to modules with at least 30 hits to ensure a large 
enough sample. Results are shown in figure [3] and 
in table [TJ compared with data from two Monte- 
Carlo studies where the tracker has been simulated 
assuming an ideal tracker geometry and after the 
alignment with data. Overall this shows that the 
alignment is already close to design specifications. 
Following the definitions of DMR, this is only an 



Figure 4: Overlap studies. The upper plot shows results in 
PXB (local u), the lower one in TIB. The modules are plot- 
ted grouped by layer and show the relative shift in the local 
u coordinate (in rcf) direction, parallel to the most sensitive 
direction, i.e. perpendicular to the strip orientation) be- 
tween overlapping module pairs. In the pixel barrel (PXB), 
the survey did not cover overlapping modules, therefore no 
improvement is visible for survey alone. In TIB, the plot 
clearly shows that survey improves the alignment. Never- 
theless, the best results were obtained after the alignment 
has been carried out. 

estimate of the modules' positions. 

3.3. Overlap studies 

There are regions of the tracker where modules 
have overlap. This reduces the effects of multiple 
scattering due to geometric reasons. The residuals 
from the two neighbouring modules, obtained in the 
same manner as in the previous method, are com- 
pared. From the results plotted in Figure [4]) it can 
clearly be seen that the alignment performs well. 

3.4. Track parameter resolution 

The previously presented results are low-level 
measures of alignment performance. To get an im- 
pression on how the tracker operates under its in- 
tended use, tracks penetrating the pixel barrel have 
been selected. Such tracks were split at the closest 
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Table 1: Results from DMR plots. RMS values of the distributions in the DMR plots (figure[3]l are given. Observe that this 
data covers the parts of the tracker hit by the cosmic ray particles. Especially in the pixel endcaps (PXE) the illumination is 
low due to the small size of the modules and the suboptimal track angles. MC simulations virere carried out using the misaligned 
and ideal geometry as starting point (column "combined" and "ideal" respectively). 
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Figure 3: Some selected plots of the DMR. The upper two plots show the distributions in the pixel barrel for the local u 
and V coordinates. Below, the plots for the TIB and TOB are shown. More details are given in table [l] 



5 



approach to the geometric center of the tracker and 
refitted as separate tracks. Then the track param- 
eters were compared at the closest approach of the 
two tracks. This procedure mimicks coUision tracks 
as if they would originate from a common vertex 
within the pixel volume. Distribution plots for all 
track parameters show that the tracker indeed per- 
forms close to design specifications. Plots for the 
distribution of px and for the impact parameters 
are shown in figure [5] 

4. Conclusions 

The studies presented here have shown that we 
are capable of aligning the inner tracker of CMS 
close to design specifications. No conclusion can 
be made for parts insufficiently illuminated by cos- 
mic rays and remaining distortion modes leaving 
invariant. Using tracks from proton collisions will 
resolve this. We are looking forward to the start of 
data taking under beam conditions, where we will 
continue our efforts to align the inner tracker as 
closely to design specifications as possible. 
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Figure 5: Resolution plots for pT and impact param- 
eters. The uppermost plot shows the distribution of the 
curvature 1/pt, the two following plots show the distribu- 
tion of the impact parameter in the xy-plane and along z re- 
spectively. All of them are compared to unaligned geometry 
and the results from a Monte-Carlo simulation. The aligned 
detector compares very well to the expected performance in 
Monte-Carlo (labelled as MC ideal). 
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